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The present application is a divisional of co-pending U.S. Patent Application 
No. 10/155,770, filed on May 25, 2002, entitled "METHODS OF AND APPARATUS 
FOR DETECTING LOW CONCENTRATIONS OF TARGET GASES IN THE FREE 
ATMOSPHERE", by Dr. Loren D. Nelson and Martin J. O'Brien (the "Parent 
Application"), priority under 35 U.S.C. 120 is hereby claimed based on the Parent 
Application, and such Parent Application is hereby incorporated by reference. 

Background of the Invention 

1. Field of the Invention 

[0001] The present invention relates generally to detecting gases, and more 
particularly to methods of and apparatus for detecting trace amounts of target gases, 
such as natural gas, remotely along long-range paths in the free atmosphere, wherein 
the exemplary natural gas detection is enabled by simultaneously remotely detecting 
methane and ethane along the same long-range path in the free atmosphere. 
2. Description of the Related Art 

[0002] In the field of gas detection, attempts are made to detect one specific, or 
"target", gas even though local conditions render such detection difficult. Such 
difficulty may be based, for example, on the fact that the target gas may only be present 
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in "trace amounts", such as one or a few parts per billion (PPB). Moreover, the target 
gas may be mixed with, or present with, water vapor and/or undesired (non-target) 
background, or "competitive", gases that may be in the same atmosphere as contains 
the target gas, for example. The background gases are referred to as "competitive" 
gases because there are overlaps in absorption spectra of the trace gases and the 
background gases, e.g., in the infrared absorption spectra of such gases. 

[0003] Many industries require facilities to detect target gases, such that there is a 
general need for accurate, fast and cost-competitive detection of target gases. 
However, the natural gas pipeline distribution system is the largest chemical 
distribution system in the United States. As a result, although equipment for detecting 
target gases other than natural gas has wide application in the United States, for 
example, the natural gas pipeline industry has the greatest need for accurate, fast, and 
cost-competitive chemical leak detection equipment and methods. This need relates in 
part to regulations that require gas utilities to perform periodic surveys for natural gas 
leaks. 

[0004] Initially, in target gas detection for the natural gas industry, there is a need to 
distinguish between natural gas as a target gas, and other combustible gases. The main 
constituents of natural gas are methane and ethane, with methane being the primary 
component. However, methane is produced by many natural biological sources, 
including animal and plant. Thus, if an elevated methane level is sensed, it does not 
necessarily mean that there is a natural gas leak. In contrast, there are no substantial 
natural ethane emission sources. However, ethane generally does not exceed twenty- 
percent of natural gas. As a result, ethane is both more difficult to detect, but is a better 
indicator of natural gas than methane. Thus, to have an optimal natural gas detector, 
there is a need for the detector to simultaneously detect both methane and ethane to 
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assure that the detected gas is from a natural gas leak and not from a natural emission 
source. 

[0005] This need to simultaneously and independently detect both ethane and 
methane is not met by current gas detection equipment. For example, flame ionization 
detectors (FID) cannot distinguish natural gas from such competitive gases. As a 
result, when currently available FID equipment is used in an attempt to detect natural 
gas, the FED equipment provides "false natural gas alarms" based on the detection of 
leaking propane tanks, leaking gasoline cans, so-called "sewer gas", and all other 
combustible gases. A natural gas pipeline utility using the FID equipment must 
respond to each false natural gas alarm although there is in fact no natural gas leak. 
Another limitation of the FID equipment is that during the detection process, it is 
generally necessary to place the equipment very close to the ground and within a 
"cloud" of the target gas that is to be detected. As a result, the tlD detection process is 
relatively slow, and FID equipment cannot be used at a place remote from the locale of 
the gas leak, for example. 

[0006] Also, Fourier Transform Infrared (FT-IR) spectro-radiometers use an 
interferometer to determine the spectral content of light passing through the free 
atmosphere. However, the output of such FT-IR instruments is based on a combination 
of all of the gases that are optically "active" (e.g., infrared absorptive) within the 
spectral region of the instrument. Thus, their temporal response is generally poor. 
Moreover, FT-IR systems are expensive, have very limited detection range through the 
free atmosphere, and cannot detect very low concentrations of target gases. 

[0007] Further, in contrast to the FID and FT-IR techniques, tunable diode laser 
absorption Spectroscopy (TDLAS), laser absorption spectroscopy (LAS), and 
differential absorption laser-based radar (DIAL) all use laser emission sources that are 
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narrow band. For example, the DIAL devices typically monitor only one or two very 
narrow spectral absorption lines. Laser-based techniques are more costly to 
manufacture, maintain and use compared to broadband techniques such as gas 
correlation radiometry (GCR), However, gas correlation radiometry (GCR) is 
generally a passive technique that relies on solar illumination or scattering, or on 
thermal emission background. Thus, GCR instruments do not have an active source of 
energy that is directed through the free atmosphere to the instrument. Further, while 
GCR instruments may be provided with filters that improve a signal to noise ratio by 
generally limiting the overall bandwidth of light admitted to a detector of the GCR 
instrument, such filters do not provide an optimized bandwidth around an optimized 
central bandpass wavelength. As a result, the sensitivity of such GCR instruments may 
be as low as 10 to 100 parts per million (PPM). 

[0008] Moreover, the FT-IR, TOLAS, LAS, DIAL and GCR technologies provide 
separate background gas and target gas channels that are interrogated sequentially. 
That is, light transmitted along a path through the free atmosphere and then through the 
background channel may be detected by a detector first. After such detection, the light 
transmitted through the same path through the free atmosphere and then through the 
target gas channel is detected by the same detector. The resulting temporal, or 
sequential, spacing of the alternating detection of the background channel and the target 
gas channel may vary from 0. 1 second to several minutes. That is, it generally takes 
more than 0. 1 seconds for these systems to provide a complete data set consisting of a 
target gas absorption measurement and an atmospheric background measurement, and 
during that time period, there may be changes in the atmospheric conditions along the 
path of the light. Thus, the light that is transmitted along the path and through the 
target gas channel may have been subjected to different atmospheric conditions along 
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the light path (e.g., atmospheric turbulence and variability) than the light transmitted 
through the background channel. As a result, the accuracy of these instruments is 
subject to a sensitivity limitation v^hen used in a dynamic atmosphere. Atmospheric 
turbulence and variability generally limit the ultimate sensitivity of these instruments in 
that the same value of instrument output provided at different times may not be based 
on the same amount of the target gas. Further, attempts to avoid such atmospheric- 
induced inaccuracies, e.g., attempts to distinguish between signals generated based on a 
target gas and on the varying atmospheric conditions, have generally been limited to 
situations in which light is transmitted only a few feet through a detection path that may 
contain the target gas to be detected. For example, it may be practical to provide 
known modulation imposed on light transmitted along a detection path that is only a 
few feet long from transmitter to detector. Given the few feet between the transmitter 
and the detector, a conductor may easily input the characteristics of the know 
modulation to the detector so that demodulation will be accurate. However, problems 
are faced in accurately demodulating the modulated light when the detection path must, 
for practical purposes, be hundreds or thousands of feet long 

[0009] In addition to the accuracy limitation due to limited sensitivity of these prior 
instruments that interrogate the separate background gas and trace gas channels 
sequentially, such instruments have limitations when attempts are made to use the 
instrument on a mobile platform, such as a truck or airborne vehicle. For example, by 
the very motion of the mobile platform, there is a dynamic, or different, atmosphere 
through which the light is transmitted. As a result, the problems of atmospheric 
turbulence and variability limit the ultimate sensitivity of the instrument used on this 
type of platform. Remote sensing systems mounted on moving platforms have the 
additional problem of variable surface reflectivity. Prior remote sensing instruments 
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mounted on mobile platforms suffered additional sensitivity limitations due to the 
temporally sequential manner of monitoring separate background and trace gas 
channels. 

[0010] Some have attempted to avoid these and other limitations of atmospheric 
turbulence and variability by isolating samples of the target gas in a closed chamber 
that is controlled to avoid turbulence and variability. However, the need to capture 
such samples from a possible location of a gas leak, and other time factors, render the 
local sampling of the target gas impractical for detecting trace amounts of target gases 
such as natural gas along miles and miles of pipeline, for example. 

[0011] Another problem faced in detecting target gases is isolation of a signal 
representing the target gas from other signals caused by background emissions. For 
example, although some types of gas detectors modulate the light just before the light 
impinges on a light receiver to remove thermal emission due to a warm instrument 
housing, such modulation at the receiver does not remove other background emissions 
such as atmospheric turbulence, jitter, beam wander, changes in the index of refraction, 
or the constant thermal emission of the atmosphere and the Earth. 

[0012] What is needed, then, in the general field of detecting gases, is a way to 
distinguish between the presence of one target gas and other gases that are normally 
present in the free atmosphere at the same time and in the same place as the target gas. 
What is also needed in target gas detection, such as for the natural gas industry, is to be 
able to distinguish between natural gas as a target gas and other combustible gases. 
Moreover, there is a need for an optimal natural gas detector that simultaneously 
detects both methane and ethane to assure that the detected gas is from a natural gas 
leak, so as to avoid fedse natural gas alarms based on the detection, for example, of the 
noted leaking propane tanks, etc. In addition, there is a need to increase detection 
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distance, that is, to increase the distance from a target gas detection instrument to a 
remote location of the target gas that is to be detected, and to also increase the detecting 
speed of such instruments. As well, there is a need to provide an instrument that will 
have a high sensitivity independently of atmospheric turbulence and variability. 
Finally, there is a need to remove undesired influences from the target gas signal so as 
to isolate detector signals representing the target gas. Such undesired influences 
include, for example, stray atmospheric fluctuations, such as humidity, atmospheric 
turbulence, changes in the index of refraction, and beam path variation; stray back light, 
such as the constant thermal emission of the atmosphere and the Earth; and system 
influences, such as jitter, beam wander and drift, and variation of source illumination. 



Summary of the Invention 

[0013] Broadly speaking, the present invention fills these needs by providing, for the 
general field of detecting gases, methods of and apparatus for distinguishing between a 
target gas and other gases that are normally in the free atmosphere at the same time and 
in the same place as the target gas. The present invention fills the needs for trace gas 
detection in the natural gas industry by an ability to distinguish between natural gas as a 
trace gas and other combustible gases. The present invention also provides a more 
optimal natural gas detector that simultaneously detects both methane and ethane to 
assure that the detected methane is from a natural gas leak, so as to avoid false natural 
gas alarms based on the detection, for example, of leaking propane tanks, etc. The 
present invention also fills these needs by substantially increasing the detection 
distance. That is, by the present invention the distance from the detection instrument to 
a location of the target gas that is to be detected may be up to about fifty-three hundred 
feet. As a result, mobile platforms, such as trucks and aircraft, may be used to carry the 
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equipment of the present invention during high-speed remote monitoring along long 
detection paths. In addition, the equipment of the present invention is provided with 
high sensitivity independently of atmospheric turbulence and variability, and the above- 
described undesired influences are removed from the trace gas signal so as to isolate 
detector signals representing the trace gas. 

[0014] An embodiment of the present invention includes a method of optimizing a 
response of a gas correlation radiometer to a trace amount of a target gas present in the 
free atmosphere along a detection path to the gas correlation radiometer. The detection 
path may also contain at least one competitive other gas the presence of which in the 
free atmosphere may interfere with detection of the trace amount of the target gas. The 
gas correlation radiometer uses an infrared filter for the response, and the method 
includes an operation of determining an absorption spectrum of the target gas modeled 
according to field parameters. Another operation determines an absorption spectrum of 
the at least one competitive gas modeled according to the field parameters. A further 
operation determines similarity and contrast between the absorption spectra of the 
atmosphere and of the target gas. Another operation determines differences between 
respective values of contrast and similarity corresponding to a plurality of band passes 
and center wavelengths of possible infra red filters. For each of the plurality of band 
passes, another operation plots the differences as a function of center wavelength. The 
infrared filter is optimized for use in the gas correlation radiometer by an operation of 
selecting a combination of infrared filter center wavelength and band pass that resuhs 
in the largest value of contrast minus similarity for the trace gas present in the free 
atmosphere along the detection path containing the at least one competitive other gas. 
A related aspect of the method involves an operation of mounting the optimized 
infrared filter in the gas correlation radiometer. 
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[0015] Another aspect of the method is that the at least one competitive gas is water 
vapor and the determining of the second absorption spectrum determines the second 
absorption spectrum corresponding to the water vapor. 

[0016] A further aspect of the method relates to the at least one competitive gas being 
water vapor and another gas, wherein the determining of the second absorption 
spectrum determines the second absorption spectrum corresponding to both the water 
vapor and the another gas. 

[0017] A yet further aspect of the method is optimizing respective responses of each 
of two gas correlation radiometers to trace amounts of the respective target gases 
ethane and methane present in the free atmosphere along the detection path to the two 
gas correlation radiometers. For the gas correlation radiometer for ethane detection the 
at least one competitive gas is a gas other than the ethane. For the gas correlation 
radiometer for methane detection the at least one competitive gas is a gas other than the 
respective methane. The method includes a further operation of performing the initial 
method once with respect to ethane as the target gas and once with respect to methane 
as the target gas so that there are provided two optimized infrared filters each having 
the selected center wavelength and bandwidth for use with the respective ethane and 
methane gas correlation radiometers to filter light transmitted through the free 
atmosphere to the respective ethane and methane gas correlation radiometers. 

[0018] A still further embodiment of the present invention includes a method of 
selecting an optimum center wavelength and an optimum bandpass of wavelengths of 
light to be processed by a gas correlation radiometer after transmission of the light 
through the free atmosphere in which there may be a trace amount of a target gas and in 
which there is likely to be at least one competitive other gas the presence of which in 
the free atmosphere may interfere with detection of the target gas. The optimum center 
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wavelength and the optimum bandpass are used in optimizing a response of a gas 
correlation radiometer to the trace amount of the target gas. The method may include 
an operation of determining a set of similarity data as a function of overlap regions 
within a spectral region. The overlap regions are for each competitive gas and the 
target gas and are those regions within the spectral region in which respective 
absorption spectra of both the target gas and the competitive gas have absorption 
characteristics. The set of similarity data include a plurality of data items within each 
of a plurality of bandpasses, wherein one of the data items corresponds to a center 
wavelength within each bandpass. The method may also include an operation of 
determining a set of contrast data as a function of non-overlap regions within the 
spectral region, the non-overlap regions being for each of the competitive gas and the 
target gas and being those regions within the spectral region in which the first 
absorption spectrum has high absorption characteristics but the second absorption 
spectrum has low absorption characteristics. The set of contrast data include a plurality 
of data items within each of a plurality of bandpasses, wherein one of the data items 
corresponds to a center wavelength within each bandpass. Optimization is completed 
by selecting the center wavelength and bandpass of an infrared filter for use with the 
gas correlation radiometer. Such selection is based on plotting a curve for each of 
various bandpasses. Each curve plots the contrast data item minus the similarity data 
item as a function of the center wavelength. The selected center wavelength and 
bandpass have the highest value of contrast data item minus the similarity data item. 

Brief Description of the Drawings 

[0019] The present invention will be readily understood by reference to the following 
detailed description in conjunction with the accompanying drawings, in which like 
reference numerals designate like structural elements. 
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[0020] Figure lA is a schematic diagram of a system of the present invention 
illustrating an apparatus for detecting trace amounts of target gases, such as natural gas, 
remotely along a fence line detection path in the free atmosphere; 

[0021] Figure IB is a schematic diagram illustrating a light source provided with a 
modulator configured with a crystal oscillator for controlling a rotating wheel that 
alternately blocks and passes a light beam from the source; 

[0022] Figures 2A and 2B depict infrared absorption spectra of a respective target gas 
and competitive gas, illustrating overlapping lines of high infrared absorption of both 
gases; 

[0023] Figures 3A and 3B combine to define a flow chart illustrating operations of a 
method of the invention for determining an optimized central wavelength and an 
optimized bandpass of an IR filter for each particular target gas; 

[0024] Figure 4 is a graph illustrating an exemplary infrared absorption spectrum for 
the target gas methane as provided by an operation of the method of Figures 3A and 
3B; 

[0025] Figure 5 is a graph illustrating a typical infrared transmission spectrum 
resulting from modeling of methane for selected parameters as provided by an 
operation of the method of Figures 3A and 3B; 

[0026] Figure 6 is a graph illustrating a typical infrared transmission spectrum 
determined for at least one competitive gas as provided by an operation of the method 
of Figures 3 A and 3B; 

[0027] Figure 7 is a graph illustrating results of an operation of the method of Figures 
3 A and 3B in which, for each of a plurality of potential IR filter bandwidths, a resultant 
value of the difference between contrast and similarity is plotted as a function of filter 
center wavelength; 
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[0028] Figure 8 is a flow chart illustrating suboperations of an operation of the 
method of Figures 3 A and 3B for determining similarity and contrast; 

[0029] Figure 9 is a graph illustrating a sub-set of a spectral region of interest in 
which three absorption curves represent atmospheric absorption in that spectral region, 
and absorption by a target gas in that spectral region, and the product of the absorptions 
of the atmosphere and methane as the target gas; 

[0030] Figure 10 is a graph illustrating a sub-set of a spectral region of interest in 
which curves represent atmospheric absorption in that spectral region, and absorption 
by a target gas in that spectral region, and the product of the atmospheric transmission 
and absorption of methane as the target gas; 

[0031] Figure 11 is a flow chart illustrating other aspects of the method of 
optimization of characteristics of an IR filter; 

[0032] Figure 12 is a schematic diagram of a system of the present invention 
configured to detect and distinguish between natural gas as a target gas and other 
combustible gases as competitive gases; 

[0033] Figure 13 is a schematic diagram illustrating lock-in amplifier circuitry of the 
system; 

[0034] Figure 14 depicts a flow chart illustrating a method of processing a target gas 
measurement technique, which relates to a ratio of measured data to null data; 

[0035] Figures 15A and 15B define suboperations of certain operations of the flow 
chart of Figure 14; 

[0036] Figure 16 is a schematic diagram illustrating circuitry for accurately 
controlling the operation of a light detector; 
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[0037] Figure 17 is a schematic diagram illustrating a vehicle-mounted system for 
detecting trace amounts of target gases, such as natural gas, remotely along a 
moderately long detection path in the free atmosphere; 

[0038] Figure 18 is a schematic diagram illustrating a stationary system for detecting 
trace amounts of target gases, such as natural gas, remotely along a long detection path 
in the free atmosphere; and 

[0039] Figure 19 is a schematic diagram illustrating an airborne system for detecting 
trace amounts of target gases, such as natural gas, remotely along a long detection path 
in the free atmosphere. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 
[0040] An invention is described for a method of and apparatus for detecting gases, 
especially trace amounts of target gases, such as natural gas, remotely along long-range 
paths in the free atmosphere. The exemplary natural gas detection is preferrably 
enabled by simultaneously remotely detecting methane and ethane along the same long- 
range path in the free atmosphere. Details are described for systems and methods in 
which an active gas correlation radiometer is configured with separate and 
simultaneously operating background and gas channels. The background channel is 
configured with a blank cell to output data that varies according to whether there is 
target gas along the path in the free atmosphere. The gas channel is configured with a 
target gas cell to output other data that is independent of whether there is target gas 
along the path in the free atmosphere. It will be obvious, however, to one skilled in the 
art, that the present invention may be practiced without some or all of these specific 
details. In other instances, well known apparatus or process operations have not been 
described in detail in order not to obscure the present invention. 
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[0041] A system 50 of the present invention is shown in Figure 1 A for detecting the 
presence of trace amounts of a specific gas 51 within the free atmosphere 52. The 
specific gas 51 is one gas of many gases that may exist in the free atmosphere 52. For 
purposes of description, the gases that may exist in the free atmosphere 52 other than 
the specific gas 51 are referred to as background, or competitive, gases 53. In the sense 
that in the system 50 the specific gas 51 is "targeted", i.e., specifically identified for 
detection, and during detection the specific gas 51 is not confused with any of the many 
other competitive gases 53, the specific gas 51 is referred to as the "target" gas. 
Similarly, in the sense that in the system 50 there may be many specific gases 51 that 
are "targeted" for detection and during detection the specific gases are not confused 
with any of the many other competitive gases 53, the specific gases 51 are referred to as 
the "target" gases. Reference herein to one target gas 51 may include reference to more 
than one target gas 51 as the context permits. The potential target gases 51 which may 
be detected by the system include, but are not limited to, CH4, NO, S02, N02, NH3, 
HN03, OH, HF, HCl, HBr, HI, CIO, OCS, H2CO, HOCl, N2, HCN, CH3C1, H202, 
C2H2, C2H6, and PH3. As described below in regard to Figure 4, such target gases 51 
have "rich" absorption spectra. 

[0042] In the absence of the system 50, local conditions in the free atmosphere 52 
may render such detection of the target gases 51 difficult. Such difficulty may be 
based, for example, on the fact that the target gases 51 may only be present in the free 
atmosphere 52 in trace amounts. Trace amounts may be in such low concentrations as 
about one or a few parts per billion (PPB), or in such concentrations as to approach a 
pure gas (-10 x 10^ PPB). Such low concentration of one or a few PPB is the lower 
end of a trace gas detection range of the system 50 with respect to certain trace gases 
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51, such as ethane, for example. Such lower end is also referred to as the "minimum 
detectable concentration" of the system 50. 

[0043] In contrast to such trace amounts, the concentration of water vapor in the free 
atmosphere 53 is generally measured as 30,000 parts-per-million (PPM)-. The water 
vapor concentration within the free atmosphere 52 is highly variable and is determined 
by the season, altitude, latitude, and local weather events, for example. Generally, such 
water vapor concentration may vary from 1,000 to 40,000 PPM, whereas the 
concentration of the target gas methane 51 is due to natural background (biological) 
sources and is typically in the range of from about 0.5 PPM to about 5.0 PPM within 
the free atmosphere 52. Since natural gas is almost entirely comprised of methane, the 
methane concentration leaking from the pipe, at the site of the leak, approaches 8 x 10^ 
PPM (-80% of pure gas). However, the methane rapidly dilutes into the free 
atmosphere 52, and the methane concentration rapidly approaches that of the natural 
background methane concentration. Due to the limited sensitivity described above, the 
Flame Ionization Detectors have a combustible gas concentration threshold on the order 
to 10 to 100 PPM to indicate a potential natural gas pipeline leak. As described below, 
the system 50 may provide a minimum detectable methane concentration of 
approximately 50 PPB, which is a higher minimum detection concentration than for 
ethane because of the lower depth of IR absorption of the methane ER absorption 
spectrum as compared to the depth of IR absorption of the ethane IR absorption 
spectrum. 

[0044] Difficulty of such detection may also result when the target gas 51 is mixed 
with, or present with, the competitive gases 53. The competitive gases 53 may include 
water vapor and/or undesired (non-target) background gases that may be in the same 
free atmosphere 52 as contains the target gas 51, for example. The competitive gases 
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53 are referred to as "competitive" gases because such gases 53 are not a target gas 51, 
yet there are overlaps 56 in absorption spectra 54 of the target gas 51 and the 
background gases 53. Figures 2 A and 2B depict respective exemplary infrared 
absorption spectra 54T and 54C of such respective target gas 51 and competitive gases 
53, wherein the overlap 56 is at an absorption line 57 having a wavenumber at about 
3015 cm V The spectrum 54C is a portion of a spectrum for United States Standard 
Atmosphere at certain conditions, and the spectrum 54T is a portion of a methane 
spectrum, for example. 

[0045] The free atmosphere 52 shown in Figure lA is an open and uncontained part 
of the Earth's atmosphere. The free atmosphere 52 is distinguished from containers 
(not shown) in which samples of a target gas 53 may be received for analysis by prior 
art equipment (not shown). Such target gas 51 in the containers is not "free". The free 
atmosphere 52 is subject to weather-related changes in pressure, temperature, 
atmospheric turbulence and variability, for example, as well as changes due to 
industrial and residential events. Such events may include discharge into the free 
atmosphere 52 of target gas 51 and competitive gases 53, for example. 

[0046] Light 58 transmitted by the system 50 may be directed along a detection path 
59 through the free atmosphere 52. At any moment of time, the light 58 is subjected to 
the current conditions of the free atmosphere 52 at a particular place, or location, along 
the detection path 59. Light 58 arriving at that same particular location at a next 
moment of time is subjected to later, then-current, conditions of the free atmosphere 52. 
Because the conditions of any such particular location of the free atmosphere 52 along 
the detection path 59 may vary from time to time (e.g., from one tenth of a second to a 
next tenth of a second, for example) the conditions of the free atmosphere 52 vary 
temporally, and are said to be dynamic. 
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[0047] Figure lA also shows the system 50 configured with a source 61 of the light 

58. The source 61 is active, that is, the source is configured to produce and direct the 
light 58 into the free atmosphere 52. The light 58 is directed along the detection path 

59, which is a beam-like region in which the target gas 51 and the competitive gases 53 
may be present. In a general sense, the light 58 is broadband radiation within a range 
of wavenumbers from about 0.2 micrometers to about one hundred micrometers. The 
source 61 of a fence line embodiment 50-1 (Figure lA) of the system 50, for example, 
may be configured to direct the light 58 from a transmission location at a near end 62 of 
the detection path 59 to a far end 63 at a receiver location spaced from the transmission 
location. Such spacing of the ends 62 and 63 may be from about thirty feet to about 
5300 feet, which in turn provides a selected length of the detection path 59 as described 
below with respect to various embodiments 50-1 through 50-5 of the system 50. 

[0048] At the transmission location, the source 61 may be provided with a modulator 
64 and beam forming optics 66, to provide the light 58 in the form of a collimated 
modulated light beam 67. The collimated modulated light beam 67 is thus generated at 
the source 61 and is transmitted from the near end 62 along the detection path 59 
through the free atmosphere 52 to the far end 63. In an embodiment 50-5 described 
below in which modulation of the light beam 67 is not required, the modulator 64 is not 
part of the source 61. The light beam 67 is transmitted along the detection path 59 
through any competitive gas 53 and through any target gas 51. During such 
transmission, depending on the atmospheric conditions along the detection path 59, the 
light 61 may be absorbed, scattered, or reflected, and such, absorption, scattering, and 
reflection may vary with location along the detection path 59, or with time, or with 
respect to the length of the detection path 59, for example. While such scattering or 
reflection may divert the light beam 67 from the detection path 59, such absorption may 
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result from infrared absorption typified by Figures 2A and 2B, for example. Such 
absorption, scattering, and reflection reduce the intensity of the light beam 67. 

[0049] Figure lA shows the light beam 67 at the far end 63 of the detection path 59 
entering a receiver 68 of the system 50. The receiver 68 is configured with a receiving 
telescope 69 that focuses the light beam 67 for transmission through an infrared (IR) 
filter 71. The IR filter 71 is configured with an optimized central wavelength and an 
optimized bandpass as described below to provide substantially increased sensitivity to 
a specific one of the target gases 51 and substantially increased selectivity of that target 
gas 51 to avoid erroneous detection of any competitive gas 53 as that target gas 51. 
The receiver 68 is configured so that the IR filter 71 transmits the light 58 having the 
now-optimized central wavelength and the now-optimized bandpass to a beam splitter 
72. The beam splitter 72 is configured to provide two split light beams 73 that are 
transmitted simultaneously to, and simultaneously enter, a dual-channel gas correlation 
radiometer (GCR) subsystem 74. 

[0050] The GCR subsystem 74 is configured with two separate channels 76 and 78 
that simultaneously receive the respective split light beams 73 having the optimized 
central wavelength and the optimized bandpass. A first of the beams 73-1 is 
transmitted through the channel 76, which is configured with a neutral density filter 79 
and a blank cell 81. The first channel 76 is configured with a first detector 82 which 
outputs blank channel data 83 that varies according to whether there is target gas 51 in 
the detection path 59 through which the light beam 67 was transmitted on the way to 
the receiver 68. A second of the beams 73-2 is transmitted through the channel 78, 
which is configured with a target gas cell 84. The second channel 78 is configured with 
a second detector 86 which outputs target gas channel data 87 that is independent of 
whether there is target gas 51 in the region, i.e. along the detection path 59 through 
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which the light beam 67 was transmitted on the way to the receiver 68. The blank 
channel data 83 and the target channel data 87 are supplied to a lock-in amplifier 88, 
the output of which is applied to a display and data processor 89. The splitting of the 
light beam 67 in this manner results in the simultaneous reception by the processor 89 
of the data 83 and 87 for simultaneous processing. 

[0051] Because of the modulation at the source 61, the modulated coUimated light 
beam 67 has been subjected to the atmospheric conditions along the detection path 59 
at successive moments of time (see successively later times tl, t2, and t3, for example, 
in Figure IB) during which the light beam 67 is transmitted from the near end 62 to the 
far end 63 of the detection path 59. Such atmospheric conditions further modulate the 
light beam 67, but this further atmospheric modulation is different from the modulation 
by the modulator 64 at the source 61. The lock-in amplifier 88 selects from this 
combination of modulated light 58 only the light 58 having the modulation imposed by 
the modulator 64 at the source 61, such that the data 83 and 87 do not include the 
effects of atmospheric turbulence, jitter, beam wander, and changes in the index of 
refraction, for example. Thus, each of the data 83 and 87 received by the processor 89 
at any moment of time represents the influence of the same current atmospheric 
conditions of the free atmosphere 52 as the light 58 is transmitted along the detection 
path 59, as may have been modified only by transmission through the respective 
channels 76 and 78. As described above, the first detector 82 outputs the blank channel 
data 83 that varies according to whether there is target gas 51 in the detection path 59 
through which the light beam 67 was transmitted on the way to the receiver 68. The 
second detector 86 outputs target gas channel data 87 that is independent of whether 
there is target gas 51 in the region, i.e. along the detection path 59 through which the 
light beam 67 was transmitted on the way to the receiver 68. The channels 76 and 78 
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thus modify the light beam 67 so that together the data 83 and 87 represent various 
factors. The factors include whether or not there was target gas 5 1 along the detection 
path 59 during the successive moments of time tl, t2, t3, etc. during which the light 
beam 67 was transmitted from the near end 62 to the far end 63 along the detection path 
59, and if target gas 51 was detected, the concentration of the target gas 51. 

[0052] Considering the structure and operation of the system 50 shown in Figure lA 
in more detail, the source 61 is configured with a broadband light source 90 that may 
output the broadband light 58 described above. The light source 90 may be a thermal 
lamp configured for high infrared output. A suitable lamp is a blackbody source having 
an output spectrum that follows the universal blackbody relationship in which the 
intensity of the light 58 output is proportional to Tcoior , the color temperature of the 
lamp. For example, the following may be used for the light source 90: silicon carbide 
lamps (having a Tcoior of -1000 Kelvin (K)), halogen lamps (having a Tcoior of -3000 
K), Xe lamps (having a Tcoior of -6000 K), tungsten lamps (having a Tcoior of -5000 K), 
and plasma glow lamps (having an output equivalent to the Xe lamp). The Xe lamp, 
for example, can achieve a total optical power output of about seventy-five Watts. The 
plasma glow lamp may be a glowing one millimeter (mm) plasma sphere, mounted at 
the focal point of an elliptical reflective mirror. The plasma glow lamp has a surface 
emission temperature of 6000 K, a total optical output of 75 watts, and an output of -6 
milliwatts within the typical bandpass of the JR filter 71. Another blackbody source 
suitable for use as the light source 90 is an arc lamp, which may be physically smaller 
and more-efficient than the glowing plasma sphere. The arc lamp produces about three 
milliwatts in the optical bandwidth (about 2800 to 3200 cm"^) required for detection of 
methane and ethane as the target gases 51. As the length of the detection path 59 
increases, it is preferable to use the source 90 having the higher T color* 
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[0053] Suitable broadband sources 61 must be capable of transmitting the generated 
light 58 into the free atmosphere 52. In order to assure that the source-generated light 
58 is transmitted into the free atmosphere 52, the broadband source 61 is not encased in 
outer bulbs, and does not have other optical elements (collectively "optical 
components"), that do not transmit well at the wavelengths of interest for the detection 
of the target gas(es) 51. When the target gas 51 is either methane or ethane, for 
example, quartz, germanium, and silicon optical components are not used since they do 
not provide for the optimal transmission of the source-generated light 58 into the free 
atmosphere 52. In a preferred embodiment of the source 61, potassium bromide may 
be used for the optical components since it has adequate optical transmission. 
However, potassium bromide scratches easily. In a more preferred embodiment of the 
source 61, sapphire is used for the optical components because sapphire transmits well 
(>90%) in a range from the visible wavelengths to approximately 7 jxm and is more 
durable than potassium bromide. 

[0054] In review, the wavelength of the light 58 output from the source 90 may be 
broadband, preferably in the range of about 0.2 micrometers to about one hundred 
micrometers, for example. More preferably, the wavelength of the light 58 output from 
the source 90 may be in the range of about 3 micrometers to about 10 micrometers. 
Most preferably, the wavelength of the light 58 output from the source 90 is in the 
range of about 3.12 micrometers to about 3.57 micrometers. Each such range of 
wavelengths may be referred to as being "broadband" in the sense that these ranges are 
substantially broader than the narrow emission spectra used in the TDLAS, laser 
absorption spectroscopy and DIAL instruments described above. 

[0055] The fence line embodiment 50-1 of the system 50 shown in Figure 1 A may be 
provided with a transmitter housing 91 at the transmitter location adjacent to the near 
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end 62 of the detection path 59. The fence line embodiment 50-1 may also be provided 
with a receiver housing 92 at the receiver location adjacent to the far end 63 of the 
detection path 59. When the detection path 59 is in a range of lengths from about thirty 
feet to about 100 feet, the modulator 64 may be hard-wired to the lock-in amplifier 88 
for directly transmitting to the amplifier 88 a demodulator signal (not shown) that is in- 
phase and in-frequency with the modulating frequency of the modulator 64. However, 
when the detection path 59 is in a range of lengths from about 100 feet to about 5,300 
feet, the modulator 64 cannot be hard-wired to the receiver 68. In this situation, as 
described below, based on receipt of data (or signals) derived from the source- 
modulated split light beams 73 the lock-in amplifier 88 determines the frequency and 
phase of the received split light beams 73 and locks to such frequency and phase. 

[0056] In one embodiment of the system 50 the source 61 may be configured with the 
modulator 64 for amplitude modulation, such as by use of a chopper having 
sequentially positioned openings, or an electro-magnetically driven tuning fork, for 
example. Figure IB shows a preferred embodiment in which the source 61 may be 
provided with the modulator 64 configured with an internal crystal oscillator 93 and 
motor control circuit 93-1 to provide closed-loop control of the modulation frequency. 
A rotating wheel 94 alternately blocks and passes the light beam 67. In this 
embodiment, the user selects the desired modulation rate through an external input 93- 
2. The motor control circuit 93-1 tunes the output of the internal oscillator 93 to this 
value. The modulator 64 includes an optical pick-off 95 associated with the wheel 94 
for providing an output signal 96-1 that is provided to the motor control circuit 93-1. 
The motor control circuit 93-1 compares the optical pick-off output 96-1 with the 
output from the crystal oscillator 93, which was previously tuned to the user input 93-2. 
Then, the motor control circuit 93-1 either speeds-up or slows-down the rotation of the 
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wheel 94 by providing appropriate drive signals 96-2 to a chopper motor 98. In this 
manner, the optical modulation frequency is precisely controlled to a user-selectable 
value using the precision crystal oscillator 93. In addition, the motor control circuit 93- 
1 provides a reference signal 97 which is precisely in phase with, and at the same 
frequency as, the optical modulation frequency. This reference signal 97 may be 
provided to the lock-in amplifier circuit 88 (Figure 13). Using this arrangement, the 
optical modulator 64 provides a precision, selectable chopping or modulation, 
frequency, which preferably is 1 kHz ±0.01Hz, or 10 parts-per-million. In a more 
general sense, the modulation frequency is well above the 1/f noise of the detectors 82 
and 86. 

[0057] In another embodiment, the source 61 may be configured with the modulator 
64 for frequency modulation, such as by use of a pressure cell through which the light 
58 from the source 90 is transmitted. 

[0058] Figure lA shows the beam forming optics 66, which may be an amateur-grade 
telescope having total reflectivity. Such telescopes are configured with totally 
reflective optics and an f-number of 4 (f/4). Such telescopes have the ability to direct 
the light beam 67 along a 5300 foot detection path 59 with a divergence of a few 
milliradians, so that if the light beam 67 is eight inches in diameter at the near end 62 
the light beam 67 may have a diameter of approximately sixty inches at the far end 63 
of the 5300 foot path 59. One advantage of the system 50 is that such telescope 66 
having the totally reflective optics for use in the IR region is very low cost as compared 
to beam forming optics constructed from refractive optics. Such optics 66 may be 
supplied by Orion, Santa Cruz, CA; or Meade of Irvine, CA, or Celestron International 
of Torrance, CA, for example. The receiving telescope 69 may also be an amateur- 
grade telescope similar to the optics 66, having the ability to receive a portion of the 
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light beam 67 having the diameter of approximately 60 inches at the end 63 of the 5300 
foot path 59. 

[0059] The IR filter 71 is configured (e.g., optimized) for each different target gas 51, 
such that the optimized central wavelength and the optimized bandpass provide 
substantially increased sensitivity to the particular target gas 51 and substantially 
increased selectivity of such target gas 51 to avoid erroneous detection of any 
competitive gas 53 as the target gas 51. In a preferred embodiment of the IR filter 71, 
the IR filter 71 is configured to respond to light wavelengths within a band 
corresponding to strong IR absorption by the specific target gas 51. The band also 
corresponds to weak IR absorption by gases 53 other than the specific target gas present 
51 in the free atmosphere 52. The word "weak" is relative to the stronger IR absorption 
by the specific target gas 51. It may be said that such band provides a high degree of 
dissimilarity in the infrared absorption spectra of the respective target gas 51 and 
competitive gas 53. 

[0060] A method of the invention is used to determine the optimized central 
wavelength and the optimized bandpass for each particular target gas 51. Figures 3 A 
and 3B show a flow chart 101 of the method. The method may start by moving to an 
operation 102 in which a modeled spectrum of the target gas 51 is determined. A 
typical initial spectrum that may be examined in operation 102 is the exemplary 
infrared (IR) absorption spectrum 103 shown in Figure 4, which is for methane, and 
extends from about wavenumber 3600 cm * to about wavenumber 700 cm'\ This 
exemplary spectrum 103 may be obtained from the spectral database "Infrared Spectra 
For Quantitative Analysis of Gases" by P. L. Hanst and S. T. Hanst, Infrared Analysis 
Inc., Anaheim, CA, for example. The target gas 51 preferably has a spectrum 103 that 
is "rich", in that there are strong absorption lines 104 in a wide range 106 of 
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wavenumbers (e.g., from about 3200 to about 1300 cm"*). The target gases 51 listed 
above are examples of target gases 51 having rich absorption spectra. 

[0061] The determination of the exemplary target gas TR spectrum 103 in operation 
102 may also include modeling the absorption features for the target gas cell 84. Such 
modeling may be performed using a suitable line-by-line database, such as the U. S. Air 
Force Research Laboratory database High Transmission Database, U.S. AI^L, 
PI7GPOS, Hanscom AFB, MA (known as "HTTRAN"). The spectral absorption 
features for a specific optical geometry may be obtained by use of computer software 
such as the U.S. AFRL FASCODE ("Fast Atmospheric Signature Code"). U.S. AFRL, 
PL/GPOS, Hanscom AFB, MA; or a PCLNWIN program made by Ontar Corporation, 
North Andover, MA, or the "HTTRANPC" program of the University of South Florida. 
Preferably, the parameters used in such modeling should generally be those to be 
experienced in the free atmosphere 52 in actual field use of the system 50. More 
preferably, the parameters should be those reasonably expected to be present in the free 
atmosphere 52 in actual field use of the system 50. Most preferably, the parameters 
used in such modeling should be those known to exist most frequently in the particular 
location of the free atmosphere 52 in which the system 50 is to be used. This modeling 
provides, or tailors, the IR absorption spectrum (e.g., the spectrum 103 of Figure 4) for 
such parameters, which may include the concentration, pressure and temperature of the 
target gas 51, and the pathlength through the optical geometry. For example, the 
wavelength, breadth and transmission characteristics of the target gas 51 change with 
changes in such parameters. Figure 5 shows a typical transmission spectrum 107 
resulting from such modeling of methane for the parameters one centimeter (cm) 
pathlength, 288.2 K methane temperature, and one atmosphere methane pressure. The 
exemplary temperature and pressure parameters are the same as the U.S. Standard 
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Atmosphere model at sea level (see Anderson, G.P, et al., Report AFGL Atmospheric 
Constituent Profiles, 0-120 km. Report No. AFGL TR-86-0110, Air Force Geophysics 
Laboratory, Hanscom AFB, MA, 1986, for example). 

[0062] The exemplary modeled absorption spectrum 107 shown in Figure 5 includes 
many deep, or highly absorbing, IR absorption lines 108D. It may be understood that 
Figure 5 represents, for methane as an exemplary target gas 51, the results of the 
operation 102, which identifies a spectral region of a first absorption spectrum of the 
target gas 51. The spectral region is the wavenumber region from 2800 to 3200 cni\ 
and the spectrum corresponds to the above-described selected parameters of target gas 
concentration, target gas temperature, target gas pressure, and path length through the 
target gas 51. The spectral region has the plurality of high (or deep) absorption lines 
108D and low absorption lines 108W. In summary, operation 102 (Figure 3A) involves 
identifying a first absorption spectrum of the target gas 51 corresponding to selected 
gas concentration, temperature, pressure, and path length. 

[0063] The method of flow chart 101 moves to operation 109 in which an absorption 
spectrum 110 (e.g., see Figure 6) is determined for at least one competitive gas 53 that 
is, and preferably for all typical competitive gases 53 that are, expected to be in the free 
atmosphere 52 in the detection path 59 in the use of the system 50. Considering the 
spectrum 110 as corresponding to all such gases 53, the spectrum 110 may be that 
obtained from a suitable line-by-line database, such as the above-identified HITRAN 
database. The spectral absorption features for a specific optical geometry of all of the 
competitive gases 53 may be obtained by use of computer software, such as one of the 
above identified FASCODE, PCLNWIN or HITRANPC programs, for example. The 
order of preference of the parameters used in such modeling are as described above for 
the target gas 51. This modeling also tailors the combined absorption spectra of the 
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competitive gases 53 for such parameters, which may include the concentration, 
pressure and temperature of the competitive gases 53, and the pathlength through the 
optical geometry. For example, the wavelength, breadth and transmission 
characteristics of the competitive gases 53 change with changes in such parameters. 
Figure 6 shows a typical absorption spectrum resulting from such modeling of the 
competitive gases 53 in the U.S. Standard Atmosphere, modeled at sea level conditions 
for the parameters of: 23 km visibility aerosols, one km pathlength, 288.2 K 
atmospheric temperature, and one atmosphere gas pressure. The exemplary modeled 
absorption spectrum 110 shown in Figure 6 also includes many deep, or highly 
absorbing, absorption lines 117D, and many weak absorbing lines 117W. The 
modeling represented by Figure 6 includes competitive gases 53 such as H20, 02, 03, 
N20, CO, C02, and CH4. Also, the U.S. Standard Atmosphere includes possible 
target gases 51 such as NO, S02, N02, NH3, HN03, OH, HF, HCl, HBr, HI, CIO, 
OCS, H2CO, HOCl, N2, HCN, CH3C1, H202, C2H2, C2H6, CH4 and PH3. 

[0064] It may be understood that Figure 6 represents, for all of the competitive gases 
53 as exemplary competitive gases 53, the results of operation 109. Figure 6 is for the 
spectral region from wavenumber 2800 to 3200 cm'\ and provides the absorption 
spectrum 110 corresponding to the above-described parameters. The spectrum 110 
includes the non-absorbing regions (the lines 117W) corresponding to the low- 
absorption characteristics (lines 108W, Figure 3A) of the exemplary methane 
absorption spectrum 107, and the highly absorbing lines 117D corresponding to deep- 
absorption characteristics (lines 108D) of such absorption spectrum 107. In summary, 
operation 109 involves identifying a second absorption spectrum of each of at least one 
other competitive gas 53. The second absorption spectrum corresponds to the selected 
gas concentration, temperature, pressure, and path length. 
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[0065] Having the absorption spectra 107 and 110, a preferred embodiment of the 
invention may compare the IR absorption spectrum 107 to the IR absorption spectrum 
110 to identify one wavelength range in which there is a high degree of dissimilarity in 
the IR absorption spectra. The dissimilarity is indicated by wavelengths in the range at 
which there is high IR absorption by the target gas 51 and low IR absorption by the at 
least one competitive gas 53. 

[0066] A more preferred embodiment of the invention configures the IR filter 71 with 
a band pass that is optimal for detecting the trace amounts of the target gas 51. The 
optimal band pass includes infrared absorption features of the target gas 51 to be 
detected and a suitable transmission region of all common gaseous constituents of the 
free atmosphere 52 other than the target gas 51. Within the band pass there is a 
relatively small degree of spectral overlap between the IR absorption features of the 
target gas IR absorption spectrum 107 and that of the atmospheric IR absorption 
spectrum 110. 

[0067] In a most preferred embodiment of the invention, the method moves to 
operation 118, which is a further part of determining optimal characteristics for 
possible IR filters 71. In a general sense, operation 118 involves determining both 
"similarity" and "contrast" that exist between the respective absorption spectra 107 and 
110 of the target gas 51 (the IR filter 71) and that of the atmosphere (i.e., the 
competitive gases 53). Once values of contrast and similarity are determined, the 
method moves to an operation 119 in which the difference between the respective 
values of contrast and similarity is determined for several center wavelengths and 
bandpasses of possible IR filters 71. 

[0068] The method moves to an operation 121 (Figure 3B) in which, for each 
potential IR filter bandpass, the resultant value of the difference between the contrast 
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and similarity is plotted as a function of filter center wavelength (see Figure 7). Each 
filter bandpass is plotted as a separate data set, that is, as a separate curve 122 in Figure 
7. Five exemplary curves 122-1 through 122-5 are shown in Figure 7, corresponding to 
five bandwidths identified respectively as 0.025 microns Half Width at Half Maximum 
C'HWHM"), 0.05 microns HWHM, 0.075 microns HWHM, 0.1 microns HWHM, and 
0.125 microns HWHM. 

[0069] The method moves to an operation 123 in which there is a selection of the 
combination of filter center wavelength and bandpass that results in the largest value of 
contrast minus similarity. Such selection provides the optimal IR filter characteristics 
for the specific trace gas 51, which is methane in the exemplary situation of Figure 5 
which shows the spectrum 107. The curves 122-1, 122-2, and 122-3, for example, 
provide choices for use in such selection. All curves 122-1 through 122-5 indicate that 
infrared filters 71 with a center wavelength from 3.2 to approximately 3.37 \im would 
not be selected due to poor performance, as indicated by the negative values of contrast 
minus similarity. That is, the target gas spectrum 107 of Figure 5 is very similar to the 
spectrum 110 of competitive gases 53 shown in Figure 6 within the spectral region 
between 3.2 to about 3.37 \im. For infrared filters 71 having a center wavelength 
between about 3.37 to 3.5 |xm, all curves 122-1 through 122-5 have positive values of 
contrast minus similarity, indicating an improved degree of contrast between the target 
gas spectrum 107 shown in Figure 5 and the spectrum 110 of competitive gases 53 
(shown in Figure 6). Within this improved region from 3.37 to 3.5 |xm, curve 122-5 
(corresponding to a filter bandpass of 0.125 fxm HWHM) provides the poorest system 
performance since it corresponds to the smallest value of contrast minus similarity. 
Alternately, curve 122-3 achieves a higher peak value of contrast minus similarity, at 
the center wavelength of approximately 3.45 \im. Also, two filter bandpasses (curves 
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122-2 and 122-3) result in equivalent values of contrast minus similarity (see HWHM 
of respective 0.05 \im and 0.075 \im). Therefore, these two infrared filters 71 
corresponding to curves 122-2 and 122-3 will provide equivalent performance for 
detecting only the target gas 51 of interest (methane in this example) while rejecting 
light 58 having wavelengths corresponding to the strong absorption bands of all 
competitive gases 53. Since the filter bandwidth of curve 122-3 is wider (0.075 \im 
HWHM), an JR filter 71 having such center wavelength (about 3.45 |im) and 
bandwidth (0.075 Jim HWHM) will transmit a larger signal amplitude of the light 67 to 
the detectors 82 and 86 than will an IR filter 71 based on curve 122-2 (0.05 pim 
HWHM). 

[0070] Having the curves 122-1 through 122-5, for example, certain additional 
considerations may be reviewed in selecting the optimal center wavelength and optimal 
bandpass, such as manufacturability of the IR filter 71. An IR filter 71 having a very 
narrow center wavelength (e.g., less than 0.5 \im) in the visible range (300 to 800 nm) 
may be manufactured more easily and at less cost than such IR filter 71 in the infrared 
or ultraviolet range. Also, consideration may be given to the breadth of the optimal 
bandpass, which may influence manufacturability and cost to manufacture, especially 
with respect to the particular center wavelength of the IR filter 71. 

[0071] More preferably, the selection of the optimal center wavelength and optimal 
bandwidth may be made according to which combination of center wavelength and 
bandpass results in the largest value of contrast minus similarity for the trace gas 51. 
For example, a center wavelength of approximately 3.45 \xm and the either the 0.05 or 
0.075 \im bandpass of curves 122-2 or 122-3, respectively, may be considered as the 
optimal center wavelength and optimal bandpass. Having made the selection of the 
optimal center wavelength and optimal bandpass, the method is DONE. 
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[0072] In view of this method, for example, when the GCR subsystem 74 is 
configured with the IR filter 71 to respond to light 67 transmitted through the detection 
path 59, the subsystem 74 (with the IR filter 71) is configured to respond to light 
wavelengths within a band in which there is a high degree of dissimilarity in the 
atmospheric infrared absorption spectra. The dissimilarity is indicated by wavelengths 
at which there is high (or strong) infrared absorption by the target gas 51 and low (or 
weak) infrared absorption by the competitive gas 53. 

[0073] Figure 8 shows in more detail suboperations of operation 118, which involves 
a suboperation 126 for determining similarity, which is defined, in general, as those 
spectral regions where both the target gas 51 (represented by the spectrum 107, Figure 
5) and the atmosphere (represented by the spectrum 110, Figure 6), possess the 
overlapping absorption features, i.e., the absorption lines 108D and 117D. 

[0074] In more detail, suboperation 126 may involve determining a set of similarity 
data as a function of overlap regions 56 (Figures 2A and 2B) within the spectral regions 
shown in Figures 2A and 2B (see absorption spectra 54C and 54T). In Figures 5 and 6 
the overlap regions may correspond to highly absorbing IR absorption lines 108D and 
1 17D for each of respective competitive gas 53 and target gas 51, and are those regions 
within the spectral region in which the respective IR absorption spectra 107 and 110 of 
both the target gas 51 and at least one other competitive gas 53 have absorption 
characteristics. The set of similarity data includes a plurality of data items (i.e., points 
on such graphs) within each of a plurality of bandpasses that include such highly 
absorbing lines 108D and 117D. One of the data items corresponds to the center 
wavelength within each bandpass. 

[0075] In still more detail, a similarity determination involves calculating the product 
of the atmospheric IR absorption (1 minus transmission T) and the target gas cell IR 
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absorption on a wavelength-by-wavelength basis. That is, a change (delta) in similarity 
is determined by: 

Equation (3) d similarity = 7^,,,, 0^-T^tmosphere)(^-T„u . 

where Tatmosphere is the wavelength-dependent optical transmission through the 

atmosphere (such as shown in Figure 6), Tceii is the wavelength-dependent optical 

transmission through the target gas channel 78 of the GCR subsystem 74 (such as that 

shown in Figure 5 for the case of methane), Tfijter is the wavelength-dependent optical 

transmission through the IR filter 71, and delta lambda is an increment of wavelength 
within the bandwidth of the IR filter 71. 

[0076] An exemplary first operation in the calculation of Equation 3 is to perform a 
product of absorptions, as follows: 

Equation (4) Product = (1 - Tatmosphere)(l " Tcell)- 

[0077] Figure 9 illustrates a sub-set of the spectral region of interest (e.g., shows 
wavenumbers 2880 through 2900 cm"^) and shows three absorption curves 131, 133, 
and 135. Curve 131 represents the atmospheric absorption in that spectral region, curve 
133 represents the absorption by the target gas 51 in that spectral region, and curve 135 
represents the product of the absorptions of curves 131 and 133 for methane as the 
target gas 51. Figure 9 shows the spectral regions where the spectrum 110 of the 
competitive gas 53 (the atmosphere in this case) and the spectrum 107 of the target gas 
51 (methane in this example) are similar, or well correlated. A high degree of 
similarity between these spectra results in a GCR that cannot distinguish between the 
competitive gas species and the target gas of interest. 
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[0078] To obtain sinnilarity, in operation 126 Equation 3 is computed as follows: 
Equation (5) similarity ^ 7,,,, a-T^mospHere)(^-Leu W , 

where X\ and X2 denote the bandwidth B of the IR filter 71. Operation 126 multiplies 
the absorption product by the filter transmission profile Tfiiter integration is 

performed over the filter bandwidth A.2 to X\, 

[0079] Still referring to Figure 8, upon conclusion of the suboperation 126 of 
operation 118, the method moves to a suboperation 137 for the determination of 
contrast. Contrast is defined, in general, as those spectral regions where the target gas 
51 (represented by the spectrum 107, Figure 5) has deep absorption features 108D, and 
the atmosphere (represented by the spectrum 110, Figure 6), possesses the weak 
absorption features, i.e., the absorption lines 117W. Thus, the determination of spectral 
contrast is based on the regions of non-overlapping absorption between the 
atmospheric spectrum 110 and the spectrum 107 of the target gas 51. 

[0080] In a general sense, operation 137 may involve determining a set of contrast 
data as a function of respective non-overlap regions 108D and 117W within the spectral 
regions of respective Figures 5 and 6. The non-overlap regions 117W and 108D are for 
each of the respective competitive gas 53 and target gas 51 and are those regions within 
the spectral region in which the target gas absorption spectrum 107 has the high 
absorption characteristics (lines 108D) but the competitive gas absorption spectrum 1 10 
has the low absorption characteristics (117W). The set of contrast data includes a data 
item corresponding to each of many possible the bandpasses and center wavelengths of 
the infrared filter 71. 

[0081] In more detail, the contrast determination of suboperation 137 involves 
calculating the product of the atmospheric transmission T and the target gas cell 
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absorption on a wavelength-by-wavelength basis. That is, a change (delta) in contrast 
is determined by: 

Equation (6) d contrast = T^,,^^ iT^tmosphere)(^''^ceu • 

[0082] An exemplary first operation in the calculation of Equation 6 is to perform 
another product of atmospheric transmission and target gas absorption per Equation (7): 

Equation (7): Product = (Tatmosphere ) (l - Tcell ); 

where: (Tatmosphere ) ^ atmospheric transmission and (l - Tceii ) = target gas 

absorption. Figure 10 illustrates a sub-set of the spectral region of interest (e.g., shows 
wavenumbers 2880 through 2900 cm'^) and shows three absorption curves 138, 139, 
and 141. Curve 138 represents the atmospheric transmission in that region, curve 139 
represents the absorption by the target gas 51 (methane) in that region, and curve 141 
represents the product of the respective transmission of curve 138 and absorption of 
curve 139 for methane as the target gas 51. Figure 10 shows the spectral regions where 
the spectrum 110 of the competitive gas 53 (the atmosphere in this case) and the 
spectrum 107 of the target gas 53 (methane in this example) are dissimilar, or not well 
correlated. A high degree of dissimilarity (or contrast) between these respective spectra 
107 and 110 results in operation of the OCR subsystem 74 that effectively distinguishes 
between the competitive gas 53 and the target gas 51 of interest. 
[0083] To obtain a value of contrast, in operation 137 Equation 6 is computed: 

Equation (8) contrast^ 

^filter 

)(i-r„„)8/i. 

[0084] Having obtained the value of contrast in suboperation 137, the method moves 
to operation 119 (Figure 3 A). 

[0085] In summary, referring to Figure 11, a method of optimization of the 
characteristics of the IR filter 71 may be described by a flow chart 146 in which an 
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operation 147 selects the series of JR filter center wavelengths for evaluation (described 
in more detail in reference to operations 102 and 109 of Figure 3A). The method 
moves to operation 148 for determining the similarity for each such IR filter center 
wavelength and bandpass selected in operations 146, as described in operation 118 
(Figure 3 A) and 126 (Figure 8). The method moves to operation 149 for determining 
the contrast for each such IR filter center wavelength and bandpass selected in 
operation 147, as described in operations 118, and 137 (Figure 8). The method moves 
to operation 151 for determining the value of contrast minus similarity for each of the 
IR filter center wavelength and bandpass options identified in operation 147, as 
described with respect to operation 119. The method moves to operation 152 in which 
the value of contrast minus similarity is plotted for each of the filter center wavelength 
and bandwidth options identified in operation 147, the plotting being as a function of 
filter center wavelength. The plotting provides a separate one of the curves 122 for 
each IR filter bandpass. The method moves to operation 153 for identifying the 
optimal center wavelength and optimal bandpass that provides the largest value of 
contrast minus similarity for the trace gas 51, as described with respect to operation 123 
(Figure 3B). 

[0086] The optimal center wavelength and optimal bandpass are used to provide the 
IR filter 71 in the system 50 to optimize the response of the GCR subsystem 74 to a 
trace amount of a target gas 71 that may be present in the free atmosphere 52 along the 
detection path 59 to the receiver 68. The detection path 59 may also contain the 
competitive gas 53 the presence of which in the free atmosphere 52 may interfere with 
detection of the trace amount of the target gas 51. Based on the IR filter 71 used with 
the GCR subsystem 74, the GCR subsystem 74 is configured to respond to the light 58 
transmitted along the detection path 59 and having a wavelength within the optimal 
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bandpass. The optimal bandpass generally corresponds to strong absorption by the at 
least one target gas 51 and corresponds to weak absorption by the gases 53 present in 
the free atmosphere 52 as compared to the strong IR absorption by the target gas 51. 

[0087] In one embodiment of the system 50, an exemplary IR filter 71 of the receiver 
68 may be provided with respect to one competitive gas 53, such as water vapor which 
is a widely occurring gas. In this case, the absorption spectrum 110 may correspond to 
water vapor. In another embodiment of the system 50, another exemplary IR filter 71 
of the receiver 68 may be provided with respect to competitive gas 53 in the form of 
water vapor and another non-target gas 53. In this case, the absorption spectrum 110 
may correspond to water vapor and the other gas 53. 

[0088] Figure 12 shows another embodiment of the system 50 which is configured to 
detect and distinguish between natural gas as a target gas 51, and other combustible 
gases as competitive gases 53. As described above, the main constituents of natural gas 
are methane and ethane. This embodiment of the system 50 may be configured as an 
optimal natural gas detector by being configured to simultaneously detect both methane 
and ethane to assure that the detected gas is from a natural gas leak and not from a 
natural emission source. Figure 12 shows the receiver 68 portion of the system 50 as a 
methane/ethane receiver 68ME provided with the receiving telescope 69 which directs 
a portion of the light beam 67 from the far end 63 to a second beam splitter 161. The 
light beam 67 is divided into two separate and simultaneous beams 162 and 163 which 
are respectively directed through two separate IR filters 71M and 7 IE. 

[0089] An embodiment of the IR filter 71M may be configured to transmit light 58 
from the light beam 67 that was transmitted along the detection path 59 and having a 
wavenumber within a band taken from the group consisting of: from about 1200 to 
about 1400 cm , from about 2800 to about 3200 cm '\ from about 2850 to about 3000 
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cm , and from about 3000 to about 3200 cm These wavenumber bands are 
preferred wavenumber bands for methane as a target gas 51. More preferrably, an 
embodiment of the JR filter 71M may be configured to transmit light 58 from the light 
beam 67 that was transmitted along the detection path 59 and having a wavenumber 
within a band taken from the group consisting of: from about 2850 to about 3000 cm , 
and from about 3000 to about 3200 cm ^ These wavenumber bands are more preferred 
wavenumber bands for methane as a target gas 51. Most preferrably, an embodiment of 
the IR filter 71M may be configured in the manner described above with respect to 
Figures 3 A and 3B through 11 with a methane-optimal central wavelength and a 
methane-optimal bandpass to provide substantially increased sensitivity to methane as 
the target gas 51 and substantially increased selectivity of methane to avoid erroneous 
detection of any competitive gas 53 as methane. The particular methane-optimal center 
wavelength and methane-optimal bandpass resulting from use of the methods described 
with respect to Figures 3 A and 3B through 11 will depend on the field conditions 
described above, for example. 

[0090] An embodiment of the IR filter 7 IE may be configured to transmit light 58 
from the light beam 67 that was transmitted along the detection path 59 and having a 
wavenumber within a band of from about 2970 to 3005 cm ■^ More preferrably, an 
embodiment of the IR filter 7 IE may be configured with an ethane-optimized central 
wavelength and an ethane-optimized bandpass in such manner described above to 
provide substantially increased sensitivity to ethane as the target gas 51 and 
substantially increased selectivity of ethane to avoid erroneous detection of any 
competitive gas 53 as ethane. Most preferrably, an embodiment of the IR filter 71M 
may be configured in the manner described above with respect to Figures 3A and 3B 
through 1 1 with an ethane-optimal central wavelength and an ethane-optimal bandpass 
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to provide substantially increased sensitivity to ethane as the target gas 51 and 
substantially increased selectivity of ethane to avoid erroneous detection of any 
competitive gas 53 as ethane. The particular ethane-optimal center wavelength and 
optimal ethane-bandpass resulting from use of the methods described with respect to 
Figures 3 A and 3B through 1 1 will depend on the field conditions described above, for 
example. 

[0091] The receiver 68ME is configured so that the IR filter 71M transmits the light 
beam 162 having, for example, the methane-optimal central wavelength and the 
methane-optimal bandpass. The light beam 162 is transmitted to a methane-ethane 
embodiment of the beam splitter 72. The beam splitter 72 is configured in two parts 
72-1 and 72-2. Each part provides two split light beams 73, including one set of split 
beams 73-1 and 73-2, and another set of split light beams 73-3 and 73-4, that are 
transmitted respectively and simultaneously to, and simultaneously enter, one of two of 
the dual channel gas correlation radiometer (GCR) subsystems 74. One subsystem 
74M for methane detection receives the beams 73-1 and 73-2; and one subsystem 74E 
for ethane detection receives the beams 73-3 and 73-4. Each respective subsystem 74M 
and 74E is configured in the manner described above with respect to Figure lA with 
two separate channels 76 and 78 that simultaneously receive the respective split light 
beams 73. Referring to Figure 1 A for the details of the respective subsystems 74M and 
74E, the first channel 76 of the methane subsystem 74M is configured with the first 
detector 82 that outputs blank channel data 83 (referred to here as 83M, for methane) 
that varies according to whether there is methane 51 in the detection path 59 through 
which the light beam 67 was transmitted on the way to the receiver 68ME. A second of 
the beams 73-2 is transmitted through the second channel 78 of the methane subsystem 
74M, which is configured with a methane target gas cell 84. The second channel 78 is 



Docket No: OPHIPOOID/CEM 



38 



PATENT 



configured with the second detector 86 which outputs target gas channel data 87 
(referred to here as 87M, for methane) that is independent of whether there is methane 
target gas 51 in the region, i.e. along the detection path 59 through which the light 
beam 67 was transmitted on the way to the receiver 68ME. The blank channel data 
83M and the target channel data 87M are supplied from the methane subsystem 74M to 
the lock-in amplifier 88, the output of which is applied to the display and data processor 
89. 

[0092] The first channel 76 of the ethane subsystem 74E is configured with a first 
detector 82 that outputs blank channel data 83 (referred to here as 83E, for ethane) that 
varies according to whether there is ethane 51 in the detection path 59 through which 
the light beam 67 was transmitted on the way to the receiver 68. The second channel 
78 of the ethane subsystem 74E is configured with a second detector 86 which outputs 
target gas channel data 87 (referred to here as 87E, for ethane) that is independent of 
whether there is ethane target gas 51 in the region, i.e. along the detection path 59 
through which the light beam 67 was transmitted on the way to the receiver 68. The 
blank channel data 83E and the target channel data 87E are supplied from the ethane 
subsystem 74E to the lock-in amplifier 88, the output of which is applied to the display 
and data processor 89. It may be understood that the splitting of the light beam 67 in 
this manner by the beam splitter 161 results in the simultaneous reception by the 
processor 89 of the data 83M, 83E, 87M and 87E for simultaneous processing to 
simultaneously detect the occurrence of both methane and ethane to assure that the 
detected gas 51 is from a natural gas leak and not from a natural emission source. 

[0093] It may be understood that the features of the embodiment 50-1 of the system 
50 may be adapted for separate, non-simultaneous detection of ethane and methane as 
the target gas 51. For example, when conditions of the free atmosphere 52 are stable 
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(relatively invariable as to atmospheric turbulence, etc.), one system 50 (as shown in 
Figure 1 A or in Figure 18 as described below) may be provided to detect ethane as a 
target gas 51 and another such system 50 may be provided to detect methane as a target 
gas 51. The detection paths 59 of each such system 50 may be close to each other (e.g., 
within a few inches or a foot). Alternatively, in the system 50 of Figure lA one dual 
channel GCR subsystem 74 may be provided to detect ethane and an additional dual 
channel GCR subsystem 74 may be combined with the ethane system 50 to detect 
methane. For example, the beam splitter 72 may be movable to direct the two beams 
73-1 and 73-2 first into the ethane subsystem 74 and then into the methane subsystem 
74. 

[0094] As described above with respect to Figures lA and IB, the modulated 
collimated light beam 67 has been subjected to the atmospheric conditions along the 
detection path 59 (e.g., at the successive moments of time tl, t2, t3, etc.) during which 
the light beam 67 is transmitted from the near end 62 to the far end 63 of the detection 
path 59. Such atmospheric conditions further modulate the light beam 67, but this 
further atmospheric modulation is different from the modulation by the modulator 64. 
The lock-in amplifier 88 selects from this combination of modulated light 58 only the 
light having the modulation imposed by the modulator 64, such that the data 83 and 87 
do not include the effects of atmospheric turbulence, jitter, beam wander, or changes in 
the index of refraction, for example. Thus, each item of the data 83 and 87 received by 
the processor 89 at any moment of time represents the influence of the same current 
atmospheric conditions of the free atmosphere 52 as the light 58 is transmitted along 
the detection path 59, as may have been modified only by transmission through the 
respective channels 76 and 78. As described above, the first detector 82 outputs the 
blank channel data 83 that varies according to whether there is target gas 51 in the 
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detection path 59 through which the light beam 67 was transmitted on the way to the 
receiver 68. The second detector 86 outputs target gas channel data 87 that is 
independent of whether there is target gas 51 in the region, i.e. along the detection path 
59 through which the light beam 67 was transmitted on the way to the receiver 68. The 
channels 76 and 78 thus modify the light beam 67 so that together the data 83 and 87 
may be processed by the processor 89 to indicate, for example, whether or not there 
was target gas 51 along the detection path 59 during the successive moments of time 
during which the light beam 67 was transmitted from the near end 62 to the far end 63 
along the detection path 59. 

[0095] The lock-in amplifier 88 is configured in two parts, one for the first channel 76 
and one for the second channel 78. Referring now to Figure 13, and considering the 
lock-in amplifier 88-1 for the first channel 76 as an example of the structure and 
operation for both channels 76 and 78, a crystal oscillator and phase lock-loop circuit 
171-1 locks to the frequency and phase of an input in the form of the blank channel 
data (or signal) 83 from the first channel 76. The crystal oscillator and lock-in circuit 
171-1 generates an output 173 modulated at the appropriate frequency. The output 173 
is delivered to a mixer 176. In addition, the lock-in amplifier 88-1 accepts a selectable 
input 172 from either the user (from knowledge of the approximate source modulation 
frequency) or from the source modulation reference signal 97 (Figure IB). A phase- 
lock-loop circuit 174-1 of the amplifier 88-1 tunes the frequency of an output 175-1 to 
that of the selectable input 172. A phase shifting circuit 174-2 within the amplifier 88- 
1 adjusts the electronic phase of the signal 175-1 from the phase lock loop circuit 174- 
1. An output 175-2 from the phase shifting circuit 174-2 is then provided to the mixer 
176, The mixer 176 within the amplifier 88-1 produces the sum and difference 
frequencies of inputs to the mixer 176 (i.e., outputs 173 and 175-2). An output 177 
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from the mixer 176 is passed through a low pass filter 178 to remove the sum of the 
frequencies contained in the signals 175-2 and 173 and thereby passing only the 
difference frequencies. Finally, the signal 177 that passes through the low pass filter 

178 is amplified by an amplifier 178-1, resulting in a lock-in output signal 179. Thus, 
the amplified signal 179 is proportional to the magnitude of the signal 83 from the 
blank channel 76 that is at the same frequency as the selectable input 172. The signal 

179 is input to the processor 89. Minor changes to the selectable input 172 frequency 
may be performed under microprocessor control in order to maximize the output signal 
179 by varying the phase lock loop 174-1 frequency via control line 180. Phase 
adjustments may also be performed under microprocessor control in order to maximize 
the output signal 179. A signal proportional to the phase adjustment performed by the 
microprocessor may be obtained at output 181. Thus, the lock-in amplifier 88-1 
outputs a signal 181 representing the phase and the signal 179 representing the 
magnitude of the blank channel signal 83 relative to the internal reference. Given this 
measured phase and magnitude of the blank channel signal 83, the lock-in amplifier 88- 
1 precisely electronically locks to the phase and frequency of the blank channel signal 
83. In a similar manner, the lock-in amplifier 88 for the first channel 78 measures the 
phase and magnitude of the target channel data 87 and precisely electronically locks to 
the phase and frequency of the target channel data 87. In this manner, (including when 
the source modulation reference signal 97, Figures IB, 17, and 18) cannot be used to 
provide the needed phase and frequency) the modulation of the light 58 at the source 61 
enables the modulated light beam 67, itself, via the blank channel data 83 which was 
derived from the light beam 67, to provide the required phase and frequency to the 
lock-in amplifier 88 of the receiver 68. In turn, this eliminates the need to extend 
cabling, or hard wire, between the modulator 64 at the source 61 and the receiver 68. 
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Further, the provision of the lock-in amplifier 88 in the receiver 68, in conjunction with 
the configuration of the dual channel GCR subsystem 74 for simultaneous channel 
operation, effectively eliminates from the first and second channel data 83 and 87 the 
effects on the light beam 67 of atmospheric turbulence, beam wander and drift, source 
illumination variations, etc. In attempts by other than the present invention to detect 
target gases 51, these atmospheric and systematic effects often dominate the target gas 
measurement uncertainty. Thus, the above-described configurations of the source 61 
with the modulator 64, and of the receiver 68 with the lock-in amplifier 88, and of the 
GCR subsystem 74, greatly improve the detectable concentration achievable by the 
system 50, which may be reduced to the above-described few PPB for ethane detection, 
for example. 

[0096] As described above, the data 83 and 87 from the respective blank and target 
gas channels 76 and 78 are input to the processor 89. The processor 89 is progranmied 
to perform operations of measurement techniques to further reduce the impact of 
competitive gases 53 on the minimum detectable concentration achievable by the 
system 50. The programmed operations of the measurement techniques also reduce the 
impact of broadband scattering, such as produced by aerosols, along the detection path 
59. These programmed operations involve measuring a value of a "null" signal under 
similar atmospheric conditions as exist during the detection and measurement of the 
concentration of the trace gas 51. 

[0097] Figure 14 shows a flow chart 181 having operations of a method of the 
invention for performing such a measurement technique. The method relates to 
Equation (1), which may be described in more detail as the following ratio: 
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Signal 



gas 



Outputs ^""''^—e. I Signal,,^^, j 

i^a^io^ii I Signal 

Equation (9): Signalf,^^,^^ J 

where Output represents the output of the processor 89. In the ratio of Equation (9), 

Signalgas is the amplitude of the output of the detector 86 (which is the target gas data 

87) measured through the gas channel 78, and Signal^iank is the amplitude of the output 
of the detector 82 (which is the blank channel data 83) measured through the blank 
channel 76. RatiOjneasured is the ratio formed by dividing Signalgas by Signalbiank- 
Rationull is determined by removing the target gas cell 84 and the blank cell 81 from 
the respective channels 78 and 76. Under this modified channel configuration, 
Signalgas-null is the amplitude of the output of the detector 86 (the target gas data 87) 

measured through the gas channel 78, and Signalbjank-null is the amplitude of the output 

of the detector 82 (the blank channel data 83) measured through the blank channel 76. 

RatiOnuU is then calculated as the ratio of these two values. Thus, RatiOnull is based on 

the differences in the optical and electronic responses of the structural elements of the 
respective blank and target gas channels 76 and 78. Such elements include optics, the 
detector 86, and typical preamplifier and post-amplifier circuits and analog-to-digital 
converters (not shown). Each of the blank channel data 83 and target gas channel data 
87 is in this manner determined relative to such null value. 

[0098] The method of the invention shown in Figure 14 may start by moving to an 
operation 182 in which calibration of the system 50 is performed to obtain the above 
Output. In detail, in operation 182 the Equation (9) Output, i.e., the (RatiOmeasured / 
RatiOouu ) -1, is determined by o[)erating the system 50 using a series of known. 
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concentration-calibrated gas standards as the target gas 51. Preferrably, such 
calibration of the system 50 may be performed once the system 50 is placed at a 
particular detection path 59 just prior to detection of an unknown target gas 51. More 
preferrably, such calibration may be performed at the factory as a step in the 
manufacture of a unit embodying the principles of the system 50. 

[0099] The method moves to an operation 183 in which a look-up table is created 
from the above Output. For each known target gas 51, the look-up table lists the 
Output of the system 50 and corresponding values of the concentration of the known 
target gas 51 in terms of the units of "concentration x pathlength" (PPM-m or PPB-m). 

[00100] The method moves to operation 184 in which the system 50 is operated 

at null status to measure Signalg^s-nuU Signalbiank-null obtain RatiOnull- 

Operation 184 of the method of flow chart 181 is shown in more detail with reference 
to Figure 15 A. Operation 184 includes a suboperation 192 in which the first and 
second channels 76 and 78 are modified to null status and a null-calibration is 
performed. Such modification is performed by removing the target gas cell 84 from the 
second channel 78 and the blank cell 81 from the first channel 76. Also, the system 50 
is placed so as to define a particular detection path 59, and the housings 91 and 92 are 
aligned so that the light beam 67 is accurately centered as it is received by the receiving 
telescope 69. Such place may be fixed, e.g., on the ground or on a mobile platform 
(e.g., on a vehicle described below). 

[00101] The method moves to suboperation 193 in which light 67 is directed 

through these modified configurations of the channels 76 and 78, to obtain a value of 

respective Signalgas-null based on the amplitude of the output of the detector 86 (the 
target gas data 87) measured through the gas channel 78, and a value of Signal^iank-null 
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based on the amplitude of the output of the detector 82 (the blank channel data 83) 
measured through the blank channel 76. The method moves to an operation 194 in 

which a value of RatiOnull is then calculated as the ratio of these two values Signalgas- 

nuU and Signal^iani^.null- The method moves to a suboperation 195 in which the system 

50 is returned to operational status with no such bypasses. The method moves to 
operation 185 (Figure 14). 

[00102] In operation 185 the system 50 is operated to detect the target gas 51 of 

unknown concentration. In preparation for this operation, the target gas cell 84 is 
configured for detecting such particular target gas 51. Also, the system 50 remains 
placed so as to define the particular detection path 59, and the housings 91 and 92 
remain aligned so that the light beam 67 is accurately centered as it is received by the 
receiving telescope 69. The system 50 is operated and the Output of the processor 89 is 
obtained. 

[00103] The method moves to an operation 186 in which this processor Output 

of the system 50 is compared to the data in the look-up table. For example, the Output 
of the system 50 is used to find the value of the concentration of the known target gas 

51 that corresponds to the Output. This value found in the look-up table is the 
concentration of the target gas 51 of unknown concentration expressed in terms of 
"concentration x pathlength" (PPM-m or PPB-m). The method is DONE. 

[00104] Details of operation 185 of Figure 14 are shown in Figure 15B. A 
suboperation 197 operates the system 50 to measure the amplitude, or value, of the 
output of the detector 86 (the target gas data 87) through the gas channel 78 to provide 

a value of Signalgas- The method moves to a suboperation 198 in which the amplitude 

of the output of the detector 82 (the blank channel data 83) measured through the blank 
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channel 76 is obtained to provide a value of Signalbjank- The method moves to a 

suboperation 199 in which a value of RatiOmeasured is obtained by computing the value 

of a ratio formed by dividing the value of Signalgas by the value of Signalbiank- The 
method moves to an operation 201 in which the values obtained in operations 194 and 
199 for the respective RatiOnull and Ratiomeasured used to compute the value of the 

ratio of Ratiomeasured to RatiOnull . i-^vi to obtain the value of the above Output. The 

method then moves to operation 186, Figure 14. 

[00105] One advantage of providing RatiOnuii in the above Equation (9) is the 

resulting removal of measurement inaccuracies associated with the elements of the 
system 50. Specifically, RatiOnuu in Equation (9) removes the effects of alignment 
variability of the source 61 with respect to the receiver 68, i.e.,transmitter-to-receiver 
alignment variability. Also removed are variabilities of the responsivities of the blank 
channel detector 82 and the gas channel detector 86, which are in terms of volts of 
output per optical watt of input. Further removed is variability of channel electronic 
gain and variability associated with the optical characteristics of the blank channel 76 
and of the gas channel 78. 

[00106] An additional benefit of measuring Rationuii in the field as provided in 

operations 192 and 193, Figure 15A ( i.e., very close to the time at which RatiOmeasured is 
determined) is that such measuring of Rationuii eliminates the effect of broad 
atmospheric absorbers such as water vapor and aerosols. Such absorbers are broad in 
the sense that they absorb throughout the entire spectral region of interest established 
by the IR filter 71. That is, when Rationuii is measured in the field, the Output 
{RatiOmeasured I RatiOnuii ) -1 is largely independent of atmospheric water vapor 
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concentration, aerosol concentration and potentially other interfering trace gas 
absorption. 

[00107] Referring now to Figure 16, an embodiment of the system 50 may be 

configured as follows. The detector 82 for the blank channel 76 is shown in Figure 16 
to illustrate the configuration of both the detectors 82 and 86, which are the same 
except for the light input. The light beam 73-1 is input to the detector 82 and the light 
beam 73-2 is input to the detector 86 as shown in Figure lA. In Figure 16 the 
exemplary detector 82 is shown configured with a respective light detector 201 having 
a response to the light beam 73-1 that is a function of the light received (i.e., the 
intensity of the light), of the detector ambient temperature, and of the magnitude of a 
bias voltage 202 input to the detector 201. A voltage supply 203 is configured to 
provide each of the detectors 82 and 86 with the common bias voltage 202. 

[00108] As part of this embodiment of the system 50, a temperature stabilizer 

206 is configured for controlling the ambient temperature of each of the detectors 82 
and 86. This temperature control is applied during all operations of the system 50 to 
minimize error due to changes in the sensitivity of the light detectors 201 to the light 58 
of the respective beams 73, as measured in amperes of output per watt of optical input. 
The change in such sensitivity of a typical light detector 201 may be approximately 4% 
per degree C, and the stabilizer 206 is effective to maintain the detector temperature 
within a tolerance not to exceed about 0.1 degree C of a desired detector temperature. 
The stabilizer 206 limits the temperature-related drift of the detectors 82 and 86 to 
about ±0.01%, for example. 

[00109] With the voltage supply 203 and the stabilizer 206 in operation as 

described, the responses of the two detectors 82 and 86 to the respective light beams 73 
are independent of the detector ambient temperature and the bias voltage 202 input to 
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the respective detector 82 and 86, and dependent on the intensity of the respective hght 
beam 73 received by the respective detector 82 and 86. 

[00110] As described above with respect to Figure lA, the source 61 of the 

system 50 may be configured to direct the light 58 from the transmission location at the 
near end 62 of the detection path 59 to the far end 63 at the receiver location spaced 
from the transmission location. In general, such spacing of the ends 62 and 63 may be 
from about ten feet to about 5300 feet, which in turn provides a similar selected length 
of the detection path 59. Thus, in the fence line embodiment 50-1 of the system 50 
shown in Figure 1 A, for example, this detection path length may be from about 10 feet 
to about 5300 feet and the separate transmission housing 91 and receiver housing 92 are 
provided at the respective spaced ends 62 and 63. 

[00111] A moderate distance embodiment 50-2 of the system 50 is shown in 

Figure 17, in which both the transmission housing 91 and receiver housing 92 are 
provided at the same near end 62 of the detection path 59. In the embodiment 50-2 the 
source 61 may be configured to direct the light beam 67 from the transmission location 
at the near end 62 of the detection path 59 along a first length 59-1 of the detection path 
59 to a light beam reflection location 211 at which a reflective surface 212 may be 
present. Such surface 212 may be a natural surface 212N such as the ground or trees, 
for example, or a non-natural surface 212A such as that of a wall of a building, for 
example. The reflective surface 212 reflects the light beam 67 out of the first length 
59-1 for return to the receiver 68 at the near end 62 along a second length 59-2 of the 
detection path 59. The spacing of the near end 62 and the reflection location 211 may 
be from about fifteen feet to about 1500 feet, which in turn provides another selected 
length of the detection path 59. In this moderate distance embodiment 50-2 of the 
system 50 the detection path length (twice the distance from the near end 62 to the 
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reflection location 211) may be as much as about 3000 feet. A typical configuration of 
the moderate length embodiment 50-2 is to mount both the transmission housing 91 and 
the receiver housing 92 on a vehicle 216, such as a truck, so that the housings 91 and 92 
are provided at the same near end 62. The direction of the detection path 59 is typically 
away from the vehicle 216, such as away from a road on which the vehicle 216 may be 
driven, so that any target gas 51 that is away from the vehicle 216 along the detection 
path 59 to or from the surface 212 may be detected. 

[00112] The embodiment 50-3 of the system 50 is shown in Figure 18, in which 
both the transmission housing 91 and receiver housing 92 are provided at the same near 
end 62 of the detection path 59. In the embodiment 50-3 the source 61 may be 
configured with a transmitting telescope 221 for long distance direction of the light 
beam 67 from the transmission location at the one end 62 of the detection path 59 along 
a second long length 59-3 of the detection path 59 to a light beam reflection location 
222 at which a remotely located reflector 223 may be present. Such reflector 223 may 
be designed to reflect the light beam 67 out of the first length 59-3 for return to the 
receiving telescope 69 of the receiver 68 at the near end 62 along a fourth length 59-4 
of the detection path 59. The spacing of the near end 62 and the reflection location 222 
may be from about 15 to about 2600 feet, which in turn provides another selected 
length of the detection path 59. In this embodiment 50-3 of the system 50 the detection 
path length (twice the distance from the near end 62 to the reflection location 222) may 
be from about thirty feet to about 5300 feet. A typical configuration of the embodiment 
50-3 is to mount both the transmission housing 91 and receiver housing 92 fixed to the 
ground at a known distance from the reflector 223, for example. With this possible 
path length of up to about a mile, the embodiment 50-3 may be referred to as a long- 
distance embodiment. 
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[00113] Another long distance embodiment 50-4 of the system 50 is shown in 

Figure 19 as combining the vehicle-mounted aspects of embodiment 50-2 and the long 
distance detection aspects of the embodiment 50-3. In embodiment 50-4 both the 
transmission housing 91 and receiver housing 92 are provided on an airborne vehicle, 
or aircraft, 230, such as an airplane (or helicopter) at the same near end 62 of the 
detection path 59. Except for the mounting of the housings 91 and 92 on the airborne 
vehicle 230 rather than on the ground, and the primary use of the ground as a reflective 
surface 231, embodiments 50-3 and 50-4 are the same. Referring to Figure 18, the 
source 61 may be configured with a transmitting telescope 221 for long distance 
direction of the light beam 67 from the transmission location at the near end 62 of the 
detection path 59. Figure 19 shows the beam 67 directed along a second long length 
59-5 of the detection path 59 to a light beam reflection location 222 which is typically 
the ground which acts as the remotely located reflector 231. Such reflector 231 reflects 
the light beam 67 out of the first length 59-5 for return to the receiver 68 at the near end 
62 (Figure 18) along a long return length 59-6 of the detection path 59, The spacing of 
the near end 62 and the reflection location 222 may be from about 0 (i.e. when the 
aircraft 230 is on the ground) to about 2500 feet, which in turn provides another 
selected length of the detection path 59. In this long distance airborne embodiment 50- 
4 of the system 50, the detection path length (twice the distance from the near end 62 to 
the reflection location 222) may be from about 0 feet to about 5000 feet, for example. 
The direction of the detection path 59 is typically directly downward away from the 
aircraft 230 toward the ground 231 so that any target gas 51 that is away from the 
aircraft 230 along the detection path 59 may be detected. 

[00114] Because of the variable conditions along the long lengths of the 

detection paths 59 of the embodiments 50-1 through 50-4, and because of the additional 



Docket No: OPfflPOOlD/CEM 



51 



PATENT 



reflection conditions which may vary widely as the aircraft 230 flies, each of the 
systems 50 for these embodiments 50-1 through 50-4 is preferably configured with the 
above-described features listed below to provide the ability to detect a low minimum 
detectable concentration of the target gas. The low minimum detectable concentration 
may vary according to the IR absorption characteristics of the target gas 51. 
Considering the above-referenced ethane and methane gases as exemplary target gases 
51, the system 50 takes advantage of the high depth of IR absorption by ethane to 
provide a lower minimum detectable concentration of ethane than is provided for 
methane. Despite the lower depth of IR absorption by methane (than by ethane), the 
system 50 still provides a minimum detectable concentration of methane (about 50 
PPB) that is orders of magnitude lower than that of the above-described prior art 
instruments. 

[00115] The listed features include the source 61 provided with the broadband 

source of light 90, and the modulator 64 at the source 61. The light source 90 may be 
the arc lamp or a configuration of a thermal emission source (e.g., a plasma glow 
lamp). The GCR subsystem 74 is provided with the beam splitter 72 so that the split 
beams 73-1 and 73-2 are processed simultaneously to achieve improved measurement 
sensitivity and accuracy. The IR filter 71 is optimized to have the optimized central 
wavelength and the optimized bandpass as described with respect to Figures 3A and 
3B, which provide substantially increased sensitivity to the particular target gas 51 and 
substantially increased selectivity of such target gas 51 to avoid erroneous detection of 
any competitive gas 53 as the target gas 51. The lock-in amplifier 88 is configured for 
generating the amplified signal 179 proportional to the magnitude of the signal 83 from 
the blank channel 76 that is at the same frequency as the selectable input 172. The 
detectors 82 and 86 are provided with the stabilizer 206 and the voltage supply 203 to 
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improve system sensitivity. Also, as described with respect to the most preferred 
embodiment of Figures 14 and 15B, the processing method of obtaining the Output 
based on Equation (9) using the (RatiOmeasured / RcitiOnuii) - i is used to remove 
measurement inaccuracies associated with the elements of the system 50, variabilities 
of the responsivities of the blank channel detector 82 and the gas channel detector 86, 
as well as the variability of channel electronic gain and variability associated with the 
optical characteristics of the blank channel 76 and of the gas channel 78. Finally, 
RatiOnuu is measured in the field (i.e., very close to the time at which RatiOmeasured is 
determined) to eliminate the effect of the broad atmospheric absorbers. 

[00116] It may be understood that other embodiments of the system 50 may 

also be configured for use in detecting target gases 51 along a short detection path 59, 
e.g., of from about 10 feet to about 100 feet in length. Also, the housings 91 and 92 
may be stationary and in the fence line configuration shown in Figure lA. The 
embodiment 50-5 is referred to as a short detection path stationary embodiment of the 
system 50. In view of the reduced length of the detection path 59 and the lack of 
adverse reflection conditions (that may vary widely, as in the mobile embodiments 50-2 
or 50-4), one may not need to include in the short detection path stationary embodiment 
50-5 of the system 50 one or more of the last mentioned features. 

[00117] In more detail, for example, in addition to the factors of atmospheric 
conditions along the detection path 59, the length of the detection path 59 is a factor 
related to the detection, determination, and measurement of the target gas 51. The 
longer the detection path 59, the more the atmospheric conditions along the detection 
path 59 may absorb, scatter, or reflect the light beam 67 that is directed along the 
detection path 59, and the more such conditions may temporally vary the light of the 
beam 67. Thus, the relatively short detection paths 59 described with respect to 
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embodiment 50-5 may be subject to different (less severe) atmospheric conditions than 
the longer detection paths 59 described with respect to Figures 18 and 19, for example. 
The atmospheric conditions to which the respective embodiments 50-3 and 50-4 of the 
system 50 of Figures 18 and 19 may be exposed render it more necessary to have the 
features of the system 50 to obtain a minimum detectable concentration of the few PPB 
(for ethane, for example). Also, surfaces 212N (Figure 19) of naturally occurring 
objects, such as the surface of the ground, or the surfaces of buildings or other things 
made by people (e.g., the reflector 222, Figure 18), may have differing light reflecting 
characteristics. Thus, the detection path situation described with respect to Figure lA 
(having no reflective surface 212 along the detection path 59) may present different 
(less severe) detection, determination, and measurement circumstances due to the 
respective reflective surfaces 223 and 231 shown in Figures 18 and 19. Similarly, the 
long detection paths 59 shown in Figures 18 and 19 may present a combination of 
circumstances based on both the reflective and the atmospheric conditions, both of 
which present different (more severe) detection, determination, and measurement 
circumstances that render it more necessary to have the features of the system 50 to 
obtain such a minimum detectable concentration of a few PPB. 

[00118] In view of the foregoing description, it may be understood that the 

present invention fills these above-described needs by providing the described methods 
of and system 50 for distinguishing between the target gas 51 and other gases 53 that 
are normally in the free atmosphere 52 at the same time and in the same place as the 
target gas 51. The present invention fills such needs for trace gas detection in the 
natural gas industry by an ability to distinguish between natural gas as a trace gas 51 
and other combustible gases 53. Also provided is a more optimal natural gas detector 
that simultaneously detects both methane and ethane to assure that the detected 
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methane is from a natural gas leak, so as to avoid false natural gas alarms based on the 
detection, for example, of leaking propane tanks, etc. Further, the method and system 
50 substantially increase the detection distance by providing the described detection 
along the long detection paths 59. Thus, by the present invention the distance from the 
system 50 to a location of the target gas 51 that is to be detected may be up to about one 
mile. As a result, mobile platforms, such as the trucks 216 and the aircraft 230, may be 
used to carry the system 50 during high-speed remote monitoring along long detection 
paths 59. In addition, the system 50 provides high sensitivity to the target gas 51 
independently of atmospheric turbulence and variability, and the above-described 
undesired influences are removed from the Output of the processor 89 so as to isolate 
the Output representing the trace gas 51. It may also be understood that in the 
subsystem 74 the channels 76 and 78 are configured so that the simultaneously 
measured data 83 and 87 vary from each other based substantially only on the presence 
of the target gas 51 in the detection path 59 at the time the broadband modulated light 
67 is directed into the free atmosphere 52 along the detection path 52. 

[00119] Although the foregoing invention has been described in some detail for 

purposes of clarity of understanding, it will be apparent that certain changes and 
modifications may be practiced within the scope of the appended claims. Accordingly, 
the present embodiments are to be considered as illustrative and not restrictive, and the 
invention is not to be limited to the details given herein, but may be modified within the 
scope and equivalents of the appended claims. 
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